A>core (bootstrap) percolation on complex networks: 
Critical phenomena and nonlocal effects 
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We develop the theory of the fc-core (bootstrap) percolation on uncorrelated random networks 
with arbitrary degree distributions. We show that the fc-core percolation is an unusual, hybrid phase 
transition with a jump emergence of the fc-core as at a first order phase transition but also with a 
critical singularity as at a continuous transition. We describe the properties of the fc-core, explain 
the meaning of the order parameter for the fc-core percolation, and reveal the origin of the specific 
critical phenomena. We demonstrate that a so-called "corona" of the fc-core plays a crucial role 
(corona is a subset of vertices in the fc-core which have exactly fc neighbors in the fc-core). It turns 
out that the fc-core percolation threshold is at the same time the percolation threshold of finite 
corona clusters. The mean separation of vertices in corona clusters plays the role of the correlation 
length and diverges at the critical point. We show that a random removal of even one vertex from 
the fc-core may result in the collapse of a vast region of the fc-core around the removed vertex. The 
mean size of this region diverges at the critical point. We find an exact mapping of the fc-core 
percolation to a model of cooperative relaxation. This model undergoes critical relaxation with a 
divergent rate at some critical moment. 

PACS numbers: 05.10.-a, 05.40.-a, 05.50.+q, 87.18.Sn 



I. INTRODUCTION 

Random damage may crucially change the structure 
and function of a network and may even completely de- 
stroy it. The description of the destruction of the com- 
plex network architectures due to damage is a challeng- 
ing direction in the multidisciplinary science of networks. 
Remarkably, much attention was attracted to the hierar- 
chical organization of various real-world networks (the 
Internet, the WWW, cellular networks, etc.) and ex- 
tracting and indexing of their highly interconnected sub- 
structures — fc-cores, cliques, and others 0,0, [J 0] 
The question is: how does random damage change and 
destroy these substructures, in particular, the fc-cores? 

The k-core of a network is its largest subgraph whose 
vertices have degree at least fc. In other words, each of 
vertices in the fc-core has at least fc nearest neighbors 
within this subgraph. The fc-core of a graph may be ob- 
tained in the following way (the fc-core algorithm or "the 
pruning rule"). Remove from the graph all vertices of 
degree less than fc. Some of the rest vertices may remain 
with less than fc edges. Then prune these vertices, and 
so on until no further pruning is possible. The result, if 
it exists, is the fc-core. The notion of the fc-core 0, 
is a natural generalization of the giant connected com- 
ponent in the ordinary percolation 0, [jj, IToL ITTl IT^L IT^ | . 
The fc-core percolation implies the breakdown of the giant 
fc-core at a threshold concentration of vertices or edges 
removed at random. In physics, the fc-core percolation 
(the fc-core percolation) on a random Bethe lattice was 
introduced in Ref. |14j for describing some magnetic ma- 
terials. The fc-core percolation on a random Bethe lattice 
was used as a basic model of the quadrupolar ordering in 
solid (o— Jjzjxjp— B^i-a mixtures ^5j], the r igid ity per- 
colation [l6j, the jamming in sphere packing [17|, glassy 



dynamics etc. An exact threshold of the emergence 
of the fc-core in some basic random networks was found 
in Ref. 0, l20| . Recently, a general exact solution of 
the fc-core problem for damaged and undamaged uncor- 
related networks with arbitrary degree distributions was 
obtained in Ref. 

These investigations revealed that the fc-core perco- 
lation is featured by an unusual phase transition which 
differs strongly from the ordinary percolation. The latter 
emerges through a continuous phase transition occurring 
at a critical concentration p c (percolation threshold) of 
the vertex occupation probability v pi H IToL fill HI fl3| . 
(A vertex is occupied with the probability p and is re- 
moved with the complementary probability Q = 1 — p.) 
At concentrations p > p c the giant connected component 
of a network occupies a finite fraction M of the total 
number of vertices TV in the network. At p — > p c , M 
tends to zero, M oc (p — p^ 13 . The standard mean-field 
exponent = 1 takes place in networks with a rapidly de- 
caying degree distribution. In networks with a scale-free 
degree distribution g -7 , exponent (3 deviates from the 
mean-field value at 2 < 7 < 4 ^^l- I n these networks, 
p c = at 7 < 3. 

The fc-core percolation for fc ^ 3 demonstrates another 
behavior. When p — > p c (k), the relative size M k of the 
giant fc-core tends to a finite value Mk(p c (k)), and at 
p < p c (k) the fc-core is absent. Note that the critical 
concentration p c (k) depends on fc. In this respect, the 
fc-core percolation looks like a first-order phase transi- 
tion with a jump of M k at the critical point. However, 
for a first-order phase transition, one would expect that 
Mfc is an analytical function of p at p > p c (k). Con- 
trary to these expectations, M k shows a singular behav- 
ior IE Q III Elf: M k {p) - M k (p c (k)) cx [p - p^k)} 1 / 2 . 
Recently, a similar phase transition was observed by nu- 
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merical simulations of the random packing of frictionless 
particles (jamming transition) [22|. In the case of the 
ordinary percolation, the critical phenomena are related 
to the divergence of the mean size of finite clusters (fi- 
nite connected components). But what is the origin of 
critical phenomena for the fc-core percolation? First im- 
portant steps in resolving this question have been made 
in Refs. |l7], |2l| where the important role of a so-called 
"corona" of the fc-core was noted. The corona is a sub- 
set of vertices in the fc-core which have exactly k nearest 
neighbors in the fc-core. 

In the present paper we develop the qualitative and 
exact theories of the fc-core percolation on complex net- 
works. Our consideration is based on an exact solution 
of the fc-core percolation on uncorrelated networks with 
arbitrary degree distributions. Specifically, we use the 
configuration model — the maximally random graphs 
with a given degree distribution |2.'"Sl |. In the large net- 
work limit, in any finite neighborhood of a vertex, these 
graphs have a tree- like local structure, i.e., — without 
loops, see, e.g., Ref. Note that in tree-like net- 

works, finite fc-cores are absent, and so we discuss only 
the giant fc-core. 

In Sec.|n]we present a qualitative picture and demon- 
strate that the critical behavior at the fc-core percolation 
threshold is related to the critical behavior of the corona 
of the fc-core. In Sec. IIIII we overview an exact solution 
describing the fc-core organization. In Sec. II VI we discuss 
the statistics of edges in a network with the fc-core and 
the meaning of the order parameter for the fc-core per- 
colation. The critical behavior of the order parameter 
is considered in Sec. Using generating functions, in 
Sec. IVII we show that the fc-core percolation threshold 
is at the same time the percolation threshold for corona 
clusters. At this point the mean size of corona clusters 
diverges. The distribution of corona clusters over sizes is 
found in Sec. I VIII Specific correlations between vertices 
in the fc-core are discussed in Sec. I Villi It is demon- 
strated that the mean intervertex distance in corona clus- 
ters plays the role of the correlation length. In Sec. IIXI 
we derive exact equations which describe the evolution 
of the degree distribution in the fc-core with increasing 
concentration of randomly removed vertices. It is demon- 
strated that a removal of even one vertex may result in a 
vast damage of the fc-core. The "diameter" of the dam- 
aged region tends to infinity at the fc-core percolation 
threshold. In Sec. El we propose an exact mapping of the 
fc-core percolation to a model of cooperative relaxation. 
This model undergoes critical relaxation with a divergent 
rate at some critical moment of the evolution. 



II. RANDOM DAMAGING THE fc-CORE 

In this section we qualitatively describe the fc-core per- 
colation in an uncorrelated random network for fc > 3. 
We assume that a vertex in the original network is oc- 
cupied with a probability p. Let M k be the probability 




FIG. 1: (a) A part of the fc=3-core with a finite corona cluster 
which consists of vertices with exactly 3 edges. This cluster 
is shown as a grey region. Corona vertices are represented 
by open dots. In a tree-like network only one corona cluster 
may connect two vertices, for example, vertices i and j on 
this figure. Removal of vertex i results in pruning all vertices 
which belong to the corona cluster. As a result, vertex j loses 
one neighbor (this is the neighboring corona vertex). The 
degree of vertex j decreases by 1. (b) In a network with loops, 
two or more corona clusters may connect together a pair of 
vertices in the fc-core. (c) In networks with nonzero clustering, 
a vertex in the fc-core may be attached to a corona cluster by 
two or more edges. In the cases (b) and (c) a removal of the 
vertex i results in pruning the corona vertices and, in turn, 
the degree of the vertex j is decreased by 2. 

that a randomly chosen vertex in a given network belongs 
to the fc-core. The fc-core in itself is a random network 
which consists of vertices with degrees at least fc. Corona 
vertices, i.e. vertices with the degree fc, are distributed 
randomly in the fc-core. They occupy a finite part of the 
fc-core. We shall show in Sec. IVII that at p > p c (k), the 
corona consists of only finite clusters, so that its relative 
size in the fc-core is sufficiently small. Any vertex in the 
fc-core may have a link with a corona vertex belonging to 
a finite corona cluster. 

Let us study the change of the fc-core size when remov- 
ing vertices at random from the original network. Let 
the occupation probability p be diminished by a value 
Ap. This corresponds to a random removal of NAp ver- 
tices. We denote the corresponding decrease in the size 
of the fc-core by NAM k = NM k (p)-NM k (p-Ap). This 
change is a quantity to be found. Firstly, there is a trivial 
contribution to NAM k due to the removal of the deleted 
vertex from the fc-core: 

NSM k = NApdM k /dp = ApNM k /p. (I) 

This can be seen from Eqs. JSJ and (JSJ). 

Secondly, after removing a vertex together with its 
edges from the fc-core we must prune all other vertices 
which will occur with degrees less than fc. In fact, the 
removal of a single vertex i from the fc-core results in the 
removal of the entire corona clusters attached to vertex 
i. Note that several corona clusters may be attached to 
a vertex with degree n > fc in the fc-core. The removal of 
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the corona clusters happens due to "the domino princi- 
ple". Indeed, after removing vertex i, its nearest neigh- 
boring corona vertex loses one link with the fc-core. This 
vertex must be pruned from the fc-core, because now it 
has only k — 1 links with the fc-core. Due to this removal, 
each of second nearest neighbors of vertex i in the corona 
clusters also loses one link with the fc-core and also must 
be pruned, and so on until all vertices in the corona clus- 
ters will be pruned one by one. This process is explained 
in Fig.^a), where a part of the fc=3-core with a corona 
cluster is represented. Let N crn be the mean total size of 
all corona clusters attached to a vertex in the fc-core [2^ . 
Then the second contribution to NAM k is N5M k N crn 
which is the number of the deleted vertices in the fc-core 
multiplied by ]V cm . 

What happens with other vertices remaining in the fc- 
core after the removal of vertex i together with corona 
clusters attached to il If there is no loops, all nearest 
neighbors of the deleted vertex i and of the pruned corona 
vertices remain in the fc-core. Their degrees decrease by 
1 since each of these vertices loses one link with the fc- 
core: n — ► ri — 1 > fc. On the other hand, in networks 
with loops, due to the pruning, a vertex may lose more 
than one connection to the fc-core. Such situations are 
represented in Figs.QIb) and^c). 

Thus, in a tree-like network, 

NAM k = NSM k + N5M k N crn . (2) 
In a differential form this equation looks as follows: 

*±*=1 + N m . (3) 
amp 

We will show in Sec. IVII that corona clusters percolate 
exactly at the fc-core percolation threshold p c (fc) and that 
iV crn diverges as [p — p c (fc)] -1 / 2 . Consequently, according 
to Eq. (0, M k demonstrates a critical singularity: 

M k (p) - M k (p c (k)) cx [p-p^k)] 1 ' 2 . (4) 

In Sec. IVII we will also show that Eq. © is exact for un- 
corrected networks with an arbitrary degree distribution 
(the configuration model) in the limit N — > 00. 

III. BASIC EQUATIONS 
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FIG. 2: (a) Schematic representation of the order parameter. 
1 — R is the probability that a given end of a randomly chosen 
edge in an undamaged network is a root of an infinite (fc— 1)- 
ary subtree. R is the probability that a given end of an edge 
is not a root of an infinite (fc— l)-ary subtree, (b) Schematic 
view of vertex configurations contributing to Mk which is the 
probability that a vertex is in the fc-core [see Eqs. @ and 
@]. A vertex in a tree-like network belongs to the fc-core if 
at least fc its nearest neighbors are the roots of infinite (fc— 1)- 
ary subtrees. The symbol V shows that the number of nearest 
neighbors which are not roots of infinite (fc— l)-ary subtrees 
is arbitrary. 

the order parameter in the fc-core problem for tree-like 
networks (see below). 

Let R be the probability that a given end of an edge 
of a network is not the root of an infinite (fc— l)-ary sub- 
tree. Then, the probability M k (n) that a vertex chosen 
at random has exactly n ^ fc neighbors in the fc-core is 
given by the equation: 

M k {n)=pY,P{q)C q n R q - n {l-R) n . (5) 

Here, P(q) is the probability that a randomly chosen 
vertex in the original undamaged network has degree q. 
C™R q - n {l - R) n is the probability that a vertex with 
degree q has q — n neighbors which are not roots of infi- 
nite (fc— l)-ary subtrees and n neighbors which are roots 
of infinite (fc— l)-ary subtrees. The combinatorial multi- 
plier C™ — m\/(m — n)\n\ gives the number of ways one 
can choose n neighbors from q neighbors. A vertex be- 
longs to the fc-core if at least fc its neighbors are roots of 
infinite (fc— l)-ary subtrees. So the probability M k that 
a vertex belongs to the fc-core is equal to 



In this section we develop an exact formalism for cal- 
culating various fc-core's characteristics. This approach 
is based on our paper [jj. 

We consider an uncorrelated network with a given 
degree distribution P(q) — the so-called configuration 
model. We assume that a vertex in the network is occu- 
pied with the probability p. In this tree- like network, the 
giant fc-core coincides with the infinite (fc— l)-ary sub- 
tree. By definition, the m-ary tree is a tree, where all 
vertices have branching at least m. The introduction of 
the (fc— l)-ary tree notion allows one to strictly define 



M k = ^M k {n). (6) 

Schematic Fig. [5] explains JSJ and JBJ. Note that for the 
ordinary percolation we must set fc = 1 in this equation. 
The number of vertices in the fc-core is equal to NM k . 

The probability R plays the role of the order parameter 
in our problem. Due to using the (fc— l)-ary trees, R is 
defined in such a way that it is independent of the second 
end of the edge — whether it belongs or does not belong 
to the fc-core. 
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(a) 




<k-2 



(b) 



FIG. 3: (a) and (b) are graphic representations of Eqs. J7( 
and JSJl, respectively. In (a) the open circle with a dashed 
boundary represents an unoccupied vertex. Other notations 
are explained in the caption to Fig. 



An end of an edge is not a root of an infinite (fc— l)-ary 
subtree if at most fc— 2 its children branches are roots of 
infinite (fc— l)-ary subtrees. This leads to the following 
exact equation for R [jj: 



k-2 



n=0 



-21 



(1-RY 



(7) 



Let us explain this equation, (i) The first term, 1— p = Q, 
is the probability that the end of the edge is unoccupied, 
(ii) C l n R l ~ n {\ — R) n is the probability that if a given 
end of the edge has i children (i.e., other edges than the 
starting edge), then exactly n of them are roots of infinite 
(k— l)-ary subtrees, (i + l)P(i + l)/^i is the probabil- 
ity that a randomly chosen edge leads to a vertex with 
branching i. Z\ — qP(q) is the mean number of the 
nearest neighbors of a vertex in the graph. Thus, in the 
square brackets, we present the probability that a given 
end of the edge has exactly n edges, which are roots of 
infinite (k — l)-ary subtrees, (iii) Finally, we take into 
account that n must be at most k — 2. In an alternative 
form, Eq. Q may be written as follows, 



R = P £ 



-k-l 



E 



(i+l)P(i+l) 



C l n R l -' l (l-R) r 



(8) 



This equation shows that a given end of an edge is a root 
of an infinite k — 1-ary tree with the probability 1 — R if it 
has at least k — 1 children (we sum over n > k — 1) which 
also are roots of an infinite (fc— l)-ary subtree. Equa- 
tions J7J) and ijSJ) are graphically represented in Fig. 
Introducing a function 

$ k{R)= E E Qff-"(i-i?r, (9) 

n— i— n 

we rewrite Eq. JJJ in a concise form: 

R=l-p + p$ k (R). (10) 



If this equation has only the trivial solution R=l, there 
is no giant k-core. The emergence of a nontrivial solution 
corresponds to the birth of the giant fc-core. The fc-core 
is described by the lowest nontrivial solution R< 1. 

The structure of the fc-core is essentially determined 
by its degree distribution Pk(n): 



(11) 



Pfc (n) is the probability to find a vertex with degree n in 
the fc-core. Note that the fc-core of an uncorrelated net- 
work in itself is an uncorrelated graph, and so it is com- 
pletely described by its degree distribution Pk(n) |27| . 
We will extensively use this circumstance. The corona oc- 
cupies a fraction Pk(k) of the fc-core. Therefore, the total 
number of vertices in the corona is equal to NMkPk{k). 
The mean degree of vertices in the fc-core is 



Zlk 



Y,Pk(n)t 



fl2) 



Comparing Eqs. (|12Jl and (JHJ, we get an important rela- 
tionship between z\ki Mk and 1 — R: 



z lk M k = Zl (l - R) 2 



(13) 



Below, in Sec. IIVI we will discuss its meaning. 

In the general case the analytical solution of Eq. (7J is 
unknown. But it can be obtained numerically. By using 
this solution, one can calculate all basic characteristics of 
the fc-core structure of an arbitrary uncorrelated network 

Some results of the numerical solution of Eq. (JJJ) for 
the Erdos-Renyi graph with z\ = 10 are represented in 
Fig. More results may be found in Ref. . This figure 
displays the dependences of the sizes of the fc-core, the 
corona and the fraction of the corona vertices in the fc- 
core on the occupation probability p for several values of 
k. One can see that far from the critical point p c (k) the 
size of the corona is small in comparison to the size of 
the fc-core. However, close to p c (k) the corona occupies 
a noticeable fraction of the fc-core. 

Let us consider a network with a scale-free degree dis- 
tribution, P(q) oc q~~' '. The case 2 < 7 ^ 3 is realized 
in most important real-world networks. With 7 in this 
range, the mean number of the nearest neighbors of a 
vertex in a network, 22, diverges if N — > 00. Solving 
analytically Eq. (|10|l . we find that the size of the fc-core 
decreases with increasing fc as follows [5j: 



M, = P 2 /(3-7) (9o/fc) (7-l)/(3-7) ! 



(14) 



where go is the minimum degree in the initial (undam- 
aged) network. Vertices which belong to the fc-core, but 
do not belong to the fc+l-core, form the fc-shell of the 
size Sk = Mk - Mk+i. Using Eq. O, at fc > 1 we find: 



S fe oc (go/fc) 2/(3 " 7) . 



(15) 
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FIG. 5: Schematic representation of the three types of edges 
in a network (large circle) with the fc-core (grey central re- 
gion): (i) edges between vertices in the fc-core (links between 
two black dots), (ii) edges between vertices which do not be- 
long to the fc-core (links between two open dots), and (iii) 
edges between vertices in the fe-core and vertices which do 
not belong to the fc-core (links between black and open dots). 



FIG. 4: Dependence of the sizes of the fc-core and the corona, 
Mk and M^(fc), respectively, on the occupation probability p 
in the Erdos-Renyi graphs with the mean degree zi = 10. 
Solid lines show Mk versus p, and dashed lines show Mk(k) 
versus p at fc = 3, 5, and 6. Notice that both Mk and Mk(k) 
have a square root singularity at the fc-core percolation thresh- 
olds, but in the curves Mk^s,e{k) this singular addition is 
practically unnoticeable. Notice the non-monotonous depen- 
dence Mk (fc) • The inset shows the fraction Pk (fc) of the corona 
vertices in the fc-core versus p. 

The asymptotic behavior given by Eqs. (|14f> and i|15|) 
agrees well with an empirical analysis of the fc-core ar- 
chitecture of the Internet on the Autonomous Systems 
level 0,0. 



IV. STATISTICS OF EDGES AND THE ORDER 
PARAMETER 

Let us consider edges in an uncorrelated network with 
the fc-core. We start with the case p = 1. There are three 
types of edges: (i) edges which connect two vertices in 
the fc-core, (ii) edges connecting two vertices which do not 
belong to the fc-core, and (iii) edges connecting together 
a vertex in the fc-core and the other one which do not 
belong to the fc-core. These types of connections in a 
network are schematically shown in Fig. |3J Let Lk, Lq 
and Lo/c be the total numbers of edges of these three types 
in the network, respectively. The sum of these numbers 
gives the total number L = Nz\/2 of edges in the initial 
network: 

L a + L k + L ok = L. (16) 

The ratios Lk/L, Lq/L and Lofe / L are probabilities that a 
randomly chosen edge in the initial network is of the type 
(i), (ii) or (iii), respectively. Because Lk = NzikMk/2, 
we can rewrite Eq. I|13f) in a form: 

^± = {l-R) 2 . (17) 



This equation has a simple explanation. It shows that the 
probability to find an edge which connects two vertices 
in the fc-core is equal to the probability that both its 
ends are roots of the (fc— l)-ary tree, that is, (1 — R) 2 
[see Fig.EJa)]. On the other hand, Eq. (|17fl explains the 
meaning of the order parameter R via the relationship 
with the measurable parameters: 1 — R = y Lk/L. 

One should note that Eq. (|17(l is also valid at p < 1 
since it follows from the exact equation ljl3|l . In this 
case, Lk must be replaced by the number Lk(p) of edges 
in a damaged fc-core while L remains the total number 
of edges in the initial network. 

Let us find the probability Lq/L that an edge chosen 
at random in the network connects two vertices which 
do not belong to the fc-core. We stress that L /L is not 
equal to R 2 as one could naively expect, but is larger, see 
Fig. Elb). Indeed, in addition to configurations where 
both the ends of an edge are not the roots of infinite 
(fc— l)-ary trees — the R 2 contribution — one must take 
into account extra configurations. In these additional 
configurations, one end of the edge is not the root of an 
infinite (fc— l)-ary tree, but the second end has exactly 
fc — 1 childrcns which are roots of infinite (fc— l)-ary trees. 
This second vertex does not belong to the fc-core as it 

(a) • • = oo oo 




(b) o O = I 1 + 2 £ | 



FIG. 6: Schematic representations of the probabilities that 
an edge connects together vertices of the fc-core or that it 
connects vertices outside the fc-core, figures (a) and (b), re- 
spectively. 
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should be. Thus we have 



L /L = R 2 + 2RY / 2 ^C q k Z 1 1 R q - k (l 



R) k ~ l . (18) 



q^k 



Zl 



Comparing the sum in Eq. (|18fl and the probability 
M k (k) given by Eq. JSJ at n = fc, we get 



L /L = R 2 + 2R 



kM k {k) 
zi{l-R)' 



(19) 



Equations l(16|) . (|17Jl and 1)19(1 establish nontrivial rela- 
tionships between independently measurable network pa- 
rameters: L, L k , Lq, Lqic, and M k (k). These relations 
may be used as a criterion of the validity of the tree 
ansatz for various networks. 

Let us now touch upon the case p < 1. After random 
removal vertices from an uncorrelated network, we again 
get an uncorrelated network. Therefore, at p < 1, we 
may still use the same formulas (|17|l . (|18|l . and (|19fl but 
with substituted characteristics of the damaged network 
— the number of edges, the mean degree, etc. |28j. 



This demonstrates that R has a jump at p = p c (k) [from 
R = R c to R = 1] as at an ordinary first-order phase tran- 
sition and a singular behavior as at a continuous phase 
transition The derivative dR/dp diverges, 



dR 
dp 

since at p 



1 - R 



«-b-p c (fc)r 1/2 , 



p[l-p& k (R)] 
Pc{k) + we have 

1- P $' fe (#)c<[p-Pc(fe)] 1/2 . 



(24) 



(25) 



This singularity suggests intriguing critical phenomena 
near the threshold of the fc-core percolation. 

In contrast, in networks with infinite z%, instead of the 
hybrid phase transition, the fc-core percolation becomes 
an infinite order phase transition [{| , similarly to the or- 
dinary percolation in this situation |9j. In this case the 
entire fc-core organization of a network is extremely ro- 
bust against random damage. 



VI. GENERATING FUNCTIONS FOR CORONA 
CLUSTERS 



V. fc-CORE PERCOLATION THRESHOLD 

When decreasing the occupation probability p, the fc- 
core decreases in size and disappears at a critical con- 
centration p c (k). According to Ref. p|, the critical con- 
centration p c (k) is determined by the following equation: 



Pc(k)K(Rc 



(20) 



Here, R c is a critical value of the order parameter R at 
the birth point of the fc-core. At p < p c {k), Eq. I|1U|) 
has only the trivial solution R = 1, and the giant fc-core 
does not exist. The derivative & k {R) = d<& k (R)/dR is 
determined by the following equation: 



p^'k(R) 



p 



E 



qP(q) 

Zl 



C q k Z\{q + l-k)R q - k {l- Rf- 



q^k 

k{k-l)P k {k)/z lk . 



(21) 



Using this equation, the condition l|20(l for the fc-core 
percolation threshold may be rewritten in the form: 



k(k-l)P k {k)/z lk = 1. 



(22) 



Let us consider the behavior of R near the phase tran- 
sition in an uncorrelated complex network with a fi- 
nite mean number z 2 of the second neighbors of a ver- 
tex, Z2 — J2 q q(q ~ 1)-P(<?)- At p near p c {k), i.e., at 
< p — p c {k) <C 1, Eq. has a nontrivial solution: 



Using the approach of Refs. 0, , we introduce the 
generating function Hi (x) of the probability that an end 
of a randomly chosen edge in the fc-core belongs to a finite 
corona cluster of a given size: 



H lk (x) = 1 - 



k -^ +x k -^[H lk{x) rK 

Z\ k Z\ k 



(26) 



Here kP k (k) / ' z\ k is the probability that an end of an edge 
chosen at random in the fc-core belongs to the corona. In 
turn, 1 — fcPfe(fc) / ' z\ k is the complementary probability 
that the end of the edge does not belong to the corona. 
We have H lk (l) = 1. 

We introduce the generating function Ho k (x) for the 
size of a corona cluster attached to a vertex in the fc-core: 



(27) 



Using this function, one can calculate the mean total size 
N CID of the corona clusters attached to a vertex randomly 
chosen in the fc-core: 



A^crn = 



dH ok (x) 



dx 



Differentiating Eqs. (|26|l and l|27|) over x gives 

fcP fc (fc) 



N m = 



l-k(k- l)P k {k)z~ k l 



(28) 



(29) 



Inserting Eqs. lj2U and (J^SJ) into Eq. we find that 
N cln diverges at the critical point, 



R c -R<x [p-pc(fc)] 



1/2 



(23) 



iVcrnOC [p - p^k)]- 1 ' 2 . 



(30) 
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This means that at p — p c (k) the corona is in its "per- 
colation transition threshold". Note however that the 
fc-core and its corona are absent at p < p c (k), so that 
there is no giant connected corona above this threshold, 
in contrast to the ordinary percolation. Equation l|22l) re- 
sembles the condition = 1 of the emergence of the 
giant connected component in the ordinary percolation. 

The exact derivation of Eq. © is based on the follow- 
ing steps. We differentiate Eq. © for M k over p: 

dM k _M k M k dR 

dp p (1 — R) dp 

Inserting Eqs. (JSTJ, (HU and (J2U into Eq. we get 

Eq. ©. 



VII. SIZE DISTRIBUTION OF CORONA 
CLUSTERS 

Let Ot crn (s) be the number of corona clusters of size 
s in the /c-core. Because the total number of vertices in 
the corona clusters is equal to NM k P k (k), we obtain the 
following condition: 



E 



sm crn {s)=NM k P k {k), 



(32) 



where s max is the size of the largest corona cluster. We 
introduce a function 



n fe (s) 



S^crn(g) 

NM k P k (k) : 



(33) 



which is the probability that a randomly chosen corona 
vertex belongs to a corona cluster of size s. The function 
LTfe(s) is related to the generating function G crn (x): 



_ 1 

s! 



dx s 



(34) 



where G CTn (x) = x[Hi k (x)] . There is a simple relation- 
ship between N crn and the mean size s crn of a corona 
cluster to which a randomly chosen corona vertex be- 
longs: 



s=l 



sTl k (s) = 1 



-JV ra 



(35) 



At s 1 the probability H k (s) has the usual asymptotic 
form 23: 



11* (a) « Cs 



-s/s* 



(36) 



Here C is a constant. Exponent a and the parameter 
s* = 1/ In I a;* I are determined by the type and the posi- 
tion x* of the singularity of the function Ho k (x) , nearest 
to x — 1. Solving Eq. (12(:i[l and inserting the obtained 
solution into Eq. I|27|l . we find that at p — > p c (k) + the 



generating functions H\ k (x) and G crn (x) have a square- 
root singularity: 

G crn (x) (xH lk (x) oc {l-x) 1 ' 2 . (37) 

In the case k = 3, Eq. (|26|l is solved exactly, 

1- {l-x/x*) 1 ' 2 



Hi k=3 (x) 



2ax 



(38) 



where a = 3P 3 (3)/zi^ =3 and x* = l/[4a(l - a)]. At the 
critical point, Eq. 122|) . we have 2a — 1 and, therefore, 
x* = 1. At p near p c (k), the parameter s* diverges, 



1/(1 -2a) 2 <xl/\p-p c (k)]. 



(39) 



The singularity l|37|l results in the standard mean-field 
exponent a — 3/2. At the critical point p = p c (k), the 
distribution function is 



U k (s) cxs- 3 / 2 . 



(40) 



It gives 9L. rn (s) cx H k (s)/s ~ s -5 / 2 . In scale- free net- 
works with a degree distribution P(q) ~ o r , this is valid 
for any 7 > 3. In contrast, in the ordinary percolation on 
scale-free uncorrelated networks, exponent a differs from 
the standard mean- field value 3/2 if 2 < 7 < 4 OS!. 

Let us estimate the size s max of the largest corona clus- 
ter. We use the condition that there is only one corona 
cluster with the size s > s max : 

OG OO 

J m cm ( S )ds = NM k P k {k) J n fc (s)«- 1 da = l. (41) 



Using asymptotics (|36|) . at p > p c (k) we obtain 



oc In N. 



(42) 



At the critical point p = p c {k), using the distribution 
function Hfc(s) oc s -3 / 2 , we get 



s oc yv 2 / 3 

5 max uc iv , 
Scrn OC V s max OC 7V 1/3 . 



(43) 

(44) 



Equation l|43() coincides with a result for the maximum 
size of a connected component at the birth point of the 
gian t connected component in classical random graphs 
[3CJ and in uncorrelated networks where three first mo- 
ments of a degree distribution converge j2^. However, 
relation l|43|l essentially differs from that for the maxi- 
mum size of a connected component if the third moment 
of the degree distribution diverges in the infinite network 
limit. 



VIII. THE CORRELATION LENGTH 

In this section we consider correlations between ver- 
tices in the /c-core with k > 3. Let us chose at random a 
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! = 1 



«0 «1 



"2 



FIG. 7: Diagrammatic representation of the mean number 
Ve(no,ni, . . . ,ni) of ways to reach a vertex which is at dis- 
tance £ from a vertex i = in the fc-core. The path goes 
through vertices m = 1,2, ...,£— 1 with degrees n m in the 
fc-core. 



vertex i in the fc-core. We aim to find the mean number 
Ve of vertices in the fc-core which may be reached from i 
following a path of a length I. In the configuration model 
the giant fc-core is a unique and simply connected sub- 
graph. Therefore, all I — 1 vertices on a path connecting 
i and j must also belong to the fc-core. Ve is given by the 
following relation: 



Ve(n ,ni, ...,ni) 



Pk(n )n Q Yl 



Pk(n m )n m (n m - 1) 



Zlk 



P k (nt)ni 



Zlk 



(45) 



A diagrammatic representation of V( is shown in Fig. \7\ 
Here n m is the degree of vertex m, where m = 0,1... £, 
on a path connecting i and j in the fc-core. For the sake 
of convenience we set i = and j = I. Here P k (no) is the 
probability that i has the degree no in the fc-core. The 
multiplier no gives the number of ways to reach vertex 
1 from i—0 following along any of its no edges. In the 
brackets, P k (n m )n m / z\ k is the probability that an edge 
outgoing from a vertex m — 1 leads to a vertex m of degree 
n rn . The multiplier n m — 1 gives the number of ways to 
leave this vertex. Finally, Pk(ne)ne/ z\ k is the probability 
that the final edge on the path leads to the destination 
vertex I of degree ne . 

Now it is easy to find the number N cin (£) of corona ver- 
tices which are at a distance £ from a randomly chosen 
vertex in the fc-core and belong to corona clusters at- 
tached to this vertex. We set ri\ = ni = ... = nt = k and 
sum over degree no of the starting vertex in Eq. (|45|1 . 
Using Eq. (|T2l gives 

Ncm(i) = ^ ^("o, ni = fc, ...ni = fc) 

nQ>k 

= kP k {k)e- {l ~ 1)/x . (46) 
Here we have introduced the correlation length 
1 1 



A = 



ln[p$' fe (i?)] ln[fc(fc - l)P k {k)/z lk ] ' 



(47) 



In accordance with Eq. (|25|l . at p — > p c {k) + the pa- 
rameter A diverges, 



(48) 



Summing iV crn (£) over £, we reproduce Eq. (|29|l for A^ crn : 



jv crn = Viv crn ffl = - fcPfc ^ 

^ 1 - exp[-l 



A]' 



(49) 



Let us determine the mean intervertex distance r crn (fc) 
between vertices in corona clusters. We use the quantity 
7-V(no,ni, . . . ,ne) and set no = ni = . . . = ng. = k. We 
find 



''crn (fc) 



Efci^MM, . . . ,fc) 
i 



l-exp[-l/A]' 



(50) 



At p close to p c (k), the correlation length A ^S> 1 and 
therefore r crn (fc) w A. 



IX. NONLOCAL EFFECTS IN THE fc-CORE 
PERCOLATION 

In Sec. ^ we have shown that a removal of a vertex 
from the fc-core leads to pruning corona clusters attached 
to the vertex. In this section we will demonstrate that 
a removal of even one vertex from the fc-core, fc > 3, 
influences degrees of vertices in a vast region of the fc-core 
around this vertex. Moreover, the size of this damaged 
region diverges at the critical point. 

Let NAp vertices be removed at random from the ini- 
tial network. As a result, the total number of vertices 
with degree n in the fc-core is changed, 

NAM k (n) = NM k (p, n) - NM k (p - Ap, n). (51) 

Let us find NAM k (n). With probability M k (n), a re- 
moved vertex may have degree n in the fc-core. Therefore, 
there is a trivial contribution to N AM k (n): 

NSM k (n) = NApdM k (n) /dp = NApM k (n) /p. (52) 

Removal of a vertex i in the fc-core may influence on the 
degree of a vertex j which is at a distance £ from i. If j 
is a nearest neighbor of i, then the degree n of vertex j 
will be decreased by 1. If £ > 1, then the probability of 
this effect is determined by the probability that j and i 
are connected by a chain of corona vertices. If vertex i is 
removed, then all vertices of a corona cluster attached to 
i also must be pruned from the fc-core due to the domino 
principle. As a result, vertex j loses one neighbor in 
the fc-core. Let V(£, n) be the mean number of vertices 
of degree n which are connected to a randomly chosen 
vertex i in the fc-core by a chain of corona vertices of 
length £. Removal at random NSM k vertices results in a 
decrease of M k (n) by a quantity 



N5M {2 \n) = NSM k ^ V{£, n). 



(53) 
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At the same time vertices with degree n+l may also lose 
one edge with the fc-core. After the pruning, they have n 
edges within the fc-core. This effect increases M k (n) by 
a quantity 

oo 

N6M (3) (n) = -NSM k V(£, n+l). (54) 
t=i 

Note that only in networks with loops, vertices in the 
fc-core may change their degree by 2 during the pruning. 
Thus, in a tree-like network there are only three contri- 
butions to NAM k (n): 

OO 

AAMfe(n) = N5M k (n)+N5M k ^Tv(l,ri) 

oo 

-NSM k ^V(l, n + l). (55) 

V(l, n) is given by the probability Vt, Eq. 

V(£,n) = Ve(n , ni=k,..., ni- X = k,n e = n) 

no>k 



nPk(n)e-W\ 



(56) 



Inserting this result into Eq. in the limit Ap — > 0, 

we get the main result of the present section: 



dM k (n) 
rflnp 

where 



M k (n)+rnM k (n) -r(n+l)M k (n+l), (57) 



1 - cxp[-l/A] 



k(k-l)M k (k) 



T,l%nM k (n) 



(58) 



The parameter r determines the mean size of a region in 
the fc-core which is damaged by a removal of one ver- 
tex chosen at random. One should stress that r de- 
pends on the entire degree distribution in the fc-core: 
r = r{M k (n)}. At p close to p c (k), we have r oc A. 
Therefore, at p — > p c (k), this size diverges. Interestingly, 
the parameter r is equal to the mean intervertex distance 
r crn in corona clusters given by expression 15U|) . 

In Eq. (|57|l . the index n can take the values n = 
k, k + 1, . . . , (7 cu t. The cutoff q cu t of the network's de- 
gree distribution P(q) depends on details of a specific 
network and its size N. 

Although we derived Eq. I|57|) by using heuristic argu- 
ments, this equation is exact for uncorrelated random 
graphs in the limit N — > oo. Equation (|57|l may be 
strictly derived by differentiating Eq. (J5J over p and using 
Eq. 

The set of Eqs. (|57|) with n from k to q cu t is a complete 
set of nonlinear equations which determine M k (n) as a 
function of p. The nonlinearity is due to the functional 
dependence of r on M k (n) . 




FIG. 8: Schematic picture of a relaxation process described 
by Eq. 16011 . An initial distribution M k (n,t = 0) over states 
with n — 3cut,9cut — l,...,fc relaxes into the final state 
{Mk(n) = 0} due to transitions of vertices from a state with 
degree n + 1 to a state with degree n. Here q cu t is the maxi- 
mum degree. 



Summing over n from k to q cut on the left and right 
hand sides of Eq. (|57|) . we obtain Eq. ©. If we know 
M k (n) for an initial network, i.e., at p = 1, then we 
can use Eq. (|57|l and find the evolution of M k (n,p) with 
decreasing p. Inserting Eqs. ijTTJ and (|12|l into (|13|l . we 
can determine the order parameter R as a function of p, 



R = l- 



<?cut 



?i 1 nM k (n,p) 



n—k 



1/2 



(59) 



Alternatively, we could find the order parameter R(p), 
solving Eq. J7J), and afterwards obtain M k (n) from 
Eq. p. 



X. MAPPING TO A COOPERATIVE 
RELAXATION MODEL 

Let us consider the k-cove percolation as an evolution- 
ary process. At time t = we have an initial uncorrelated 
network with the fc-core. During a time interval At we 
remove at random a fraction Ap/p — At of occupied ver- 
tices from the network. This means that the occupation 
probability p decreases in time as p — e . With this 
substitution, Eq. (|57|l takes a form 



dM k (n,t) 
dt 



-M k (n, t)-rnM k (n, t)+r(n+l)M k (n- 



•M). 
. ( 60) 

This rate equation describes the relaxation of an initial 
distribution {M k (n,t = 0)} to the final state with the 
destroyed fc-core, i.e., {M k (n) = 0}, due to the chain of 
transitions of vertices from states of degree n+l to states 
of degree lower by one: n+l —> n, see Fig. |SJ Note that 
we consider only relaxation in states with n > k, assum- 
ing that vertices of degree less than k are pruned instantly 
[31| . The parameter r plays the role of the characteristic 
scale of the relaxation rate. This relaxation is a cooper- 
ative process due to the functional dependence of r on 
M k (n,t), see Eq. l(5H)l. At time t c (k) = ln[l/p c (fc)] this 
model undergoes a dynamical phase transition. Using 
Eq. (0J , we find that the total number M k (t) of vertices 
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in the fc-core has a singular time dependence near t c (k): 

M k (t) - M k (t c (k)) <x \p(t)-p c (k)] u cx [t c (k) - t] v . (61) 

The critical exponent v = 1/2 is valid for fc 3. Inser- 
ting Eq. (|48|) into Eq. I]58fl. we find that the relaxation 
rate diverges at the critical time i c (fc), 

r cx [p(t) -p c (fc)]-V2 oc [t c (k) - t]- 1 ' 2 . (62) 

Note that in accordance with the results obtained in 
Sec. |n] the characteristic scale r of the relaxation rate 
also determines the mean size of the region in the fc-core 
cropped out due to the deletion of a vertex. In its turn, 
r is approximately equal to the correlation length A, i.e., 
the larger is the correlation length the larger is the re- 
laxation rate. This is in contrast to the usual critical 
slowing down of the order parameter relaxation for con- 
tinuous phase transitions. In the latter case, the larger is 
the correlation length the smaller is the relaxation rate, 
r ~ A _z , where z is a dynamical critical exponent. 

XI. CONCLUSIONS 

In this paper we have explained the nature of the fc- 
core percolation transition in uncorrelated networks. To 
obtain our results, we used heuristic arguments and de- 
veloped an exact theory. Let us list the main features 
of the quite unusual fc-core percolation transition: (i) a 
jump emergence of the fc-core, (ii) the critical singular- 
ity in the phase with the fc-core, (iii) the absence of any 
critical effects — strong correlations, divergent "suscep- 
tibilities" , etc. — on the "normal phase" side. We had to 
reveal the meaning of the order parameter in this prob- 
lem, to explain the nature of the jump, to find the origin 
of the singularity, to indicate a "physical" quantity, which 
diverges at the critical point, to indicate long-range cor- 
relations, specific for this transition. 

We have shown that the order parameter in this prob- 
lem is simply expressed in terms of the relative number of 
edges in the fc-core, see relation (I16f) . The tree ansatz has 
allowed us to find the fc-core order parameter and other 
fc-core characteristics of various uncorrelated networks. 

We have found that the unique properties of the fc-core 
percolation transition are essentially determined by the 
corona subset of the fc-core, that is, by vertices with ex- 
actly fc connections to the fc-core. These are the "weak- 
est" vertices in the fc-core. The critical correlations in 



the fc-core are due to the correlations in the system of 
the corona clusters. 

In the "fc-core phase", the corona clusters are finite, 
but their sizes and long-range correlations grow as the 
network approaches the fc-core percolation threshold. 
The mean size of a corona cluster to which a randomly 
chosen corona vertex belong diverges at the fc-core perco- 
lation threshold. This quantity plays the role of a critical 
susceptibility in this problem. So, the fc-core percola- 
tion threshold coincides with the percolation threshold 
for corona clusters, and the fc-core phase is the "normal" 
phase for the corona. The dramatic difference from the 
ordinary percolation is that the corona disappears on the 
other side of the threshold, and so critical fluctuations in 
the phase without the fc-core are absent. 

For understanding the nature of this transition, wc 
have studied the process of the destruction of the fc-core 
due to the random deletion of vertices. The deletion of 
a vertex in the fc-core results in the clipping out the en- 
tire adjacent corona clusters from the fc-core due to the 
domino principle. This effect is enormously increased 
when corona clusters become large — near the fc-core 
percolation threshold. In the threshold, the removal of a 
tiny fraction of vertices results in the complete collapse 
of the corona and the fc-core. In this respect, the fc-core 
percolation problem can be mapped to a model of co- 
operative relaxation, which undergoes critical relaxation 
with a divergent rate at some critical moment. 

To conclude, let us indicate a possible application — 
a social network model, where social links connect indi- 
viduals. Each of vertices (individuals) in our model may 
occur in one of a few states — distinct beliefs, opinions, 
religions, ideologies, diseases, etc. We assume that each 
vertex takes a specific state if at least fc its neighbors are 
in this state. Is it possible that in this social net a giant, 
say, religious group will emerge? The answer is yes if the 
network has the giant fc-core. A giant community of in- 
dividuals being in the same state forms the fc-core of this 
network. We believe that our results are applicable to a 
variety of complex cooperative systems of this kind. 
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